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Abstract. The Lorentz covariant theory of precise Doppler measurements (PDM) based 
on the retarded Lienard- Wiechert solution of the Einstein equations is described. An exact 
solution of equations of light propagation in the field of arbitrary moving bodies, which 
drastically extends the range of applicability of the new theory of PDM, is obtained. An 
explicit formula for the gravitational shift of frequency is given in analytic form. The limiting 
cases of the Doppler observations in gravitational lensing and of the spacecraft's Doppler 
tracking are described in more detail. We also present the post-Newtonian theory of the 
PDM developed for searching relativistic effects in close optical binaries and massive planetary 
systems. 
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1 Background 

Detection of low-frequency gravitational waves and verification of General Relativity 
in the Solar system Doppler tracking experiments, as well as searching for anysotropies 
in the cosmic background radiation (CMB) appeal to the development of the relativis- 
tic theory of precise Doppler measurements (PDM) . Previous attempts to create such 
a theory in framework of General Relativity (GR) were based on making use of the 
post-Newtonian approximations (PNA) which operate with the instantaneous values 
of metric tensor taken on hypersurfaces of constant time. Hence, an adequate treat- 
ment of propagation of light rays in the gravitational field of an isolated astronomical 
system can be achieved using the PNA only for light-time intervals which do not ex- 
ceed the characteristic Keplerian time of the system. If light propagates longer, the 
residual terms of the PNA theory of PDM become unmanageable. Moreover, the ques- 
tion about the instant of time at which positions and velocities of gravitating bodies 
deflecting light rays should be fixed to minimize calculational errors cannot be solved 
in the framework of the PNA unambiguously. 

A natural way to overcome incompleteness in the PNA approach to the theory 
of PDM is to employ the post-Minkowskian approximation (PMA) scheme, which is 
Lorentz covariant and does not involve the slow-motion expansions with respect to 
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powers of v/c, where v is the characteristic velocity of the light-ray-defiecting body, 
and c is the speed of light. We have succeeded in developing of such a Lorentz covari- 
ant theory of PDM which is based on the retarded Lienard- Wiechert solution of the 
Einstein equations, linearized with respect to Newton's gravitational constant G [Q. 
This makes it possible to obtain a unique exact solution of the equations of light propa- 
gation in the field of arbitrary moving spherically-symmetric bodies, which drastically 
extends the range of applicability of the new theory of PDM. 



2 Doppler effect in Special and General Relativity 

Special relativistic treatment of the frequency shift is well-known. It is based on two 
facts: 1) the proper time runs differently for identical clocks moving with different 
velocities; -2) electromagnetic waves propagate along straight lines in the flat space- 
time (see the paper for more details). 

General relativistic formulation of the frequency shift in curved space-time is more 
involved. Two definitions of the Doppler shift are used || - in terms of energy (A) 
and in terms of frequency (B) 

(A) - = ^~, (B) f = §, (1) 

where Uq and v are emitted and observed electromagnetic frequencies of light; T , Uq 
and T, u a are the proper time and 4-velocity of the source of light and observer; /Co a 
and K a are null 4- vectors of the light particle at the points of emission and observation 
respectively. Despite the obvious difference in the two definitions they are identical 
as u a = dx a /dT, IC a = dip/dx a , and the phase of electromagnetic wave ip is constant 
along the light rays. In GR there does not exist a parallel transport of 4-vectors in 
the sense of the flat space-time. Instead of that one has to integrate equations of light 
propagation to connect physical quantities at the points of emission and observation 
of light. 

The integration of null geodesic equations in the case of space-times possesing sym- 
metries was known for a long time and extensivly used in astronomical practice [Q. 
However, the most interesting astronomical phenomena in the propagation of light 
rays in curved space-time are caused by small, time-dependent perturbations of the 
background geometry. Usually, the 1st PNA in the relativistic N-body problem with 
fixed or uniformly moving bodies was applied to consider effects of such perturbations 
|y . Unfortunately, this approximation works properly if, and only if, time of propa- 
gation of light is much shorter than the characteristic Keplerian time of the N-body 
problem. An adequate treatment of the effects in propagation of light must account 
for the retardation effects in the propagation of gravitational field to the point of the 
field's interaction with the electromagnetic wave. 



3 Predictive Relativistic Mechanics of Photons 



Physically consistent description of gravitational field in the problem of propagation 
of light rays is achieved in the framework of the post-Minkowskian approximations 
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(PMA) of the Einstein gravitational field equations. For our purpose one PMA is 
completely enough. This approximation is linear with respect to the universal gravity 
constant G and gives metric tensor as an algebraic sum of the Minkowski metric rj a p 
and its perturbation h a p written in terms of the Lienard- Wiechert tensor potential. 
The perturbation h a /3[t,x, x a (s), v a (s)] is a function of a field point (i,x), masses m a 
as well as coordinates x (s) and velocities v a (s) of gravitating bodies taken at the 
retarded time s = t — i|x — x a (s)|. Substitution of h a p into the equations of the 
light ray geodesies bring them to the form of to the "retarded-functional differential 
system" (hereafter G = c = 1) 

N 

x(f) - maFa M*). *(*). *«(*). v "( s )] . ( 2 ) 

a=l 

because of the dependence of the gravitational light-deflecting force F a on the retarded 
time argument s. Studying of such equations belongs to the framework of "predictive 
relativistic mechanics" ||. 



4 Solution of the light-ray geodesies 

The system (|J) is transformed to a simpler form by making use of a specific differential 
identity |lj applied to the force F a in order to change of the order of operations of 
taking partial derivatives and substitution for the unperturbed light- ray trajectory. 
Solution of the system (^) is, then, achieved by performing integrals from F a along 
the light-ray path by making use of transformation from the coordinate time t to the 
retarded time s which eliminates the time t from all of the integrands of any integral. 
Hence, the integrals can be directly performed as soon as the motion of the light-ray- 
dcflccting bodies is known. As a result, the time of light travelling from the point of 
emission xo to the point of observation x is obtained in the following form Q 

A f s [l-k.v a (Q] 2 dX 

- y v^m **+k. Xa (o-c 

where t* is a pre-defined constant. General formula (|J) is used for calculation of the 
Doppler shift in what follows. 



£-£ = |x-x |+ 2^m a / L - ^ - , z J 7 , (3) 

a=l 



5 Doppler Effect in Gravitational Lenses 

In the case of a gravitational lens having a total mass M., moving with velocity V 1 with 
respect to the chosen coordinate system, and possesing a time-dependent quadrupolc 
moment 2^ the gravitational shift of frequency reads 

= 4 (-MV* &i + ) ^ |£| + ^ (v • a) + - (v • a) , (4) 

where r, tq are distances from the lens to observer and source of light, R ~ r + r*o, 
a is the angle of the total deflection of light rays, v and Vq are velocities of observer 
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and the source of light respectively, £ = £ 8 is the impact parameter of the light ray, 
di = d/dC- 



6 Spacecraft Doppler Tracking 



Lorentz covariant approach allows to derive general relativistic corrections to the well- 
known special relativistic expression for the electromagnetic frequency shift in the 
spacecraft's Doppler tracking. The corrections are summarized in the formula 



v 1 — k • v 
v 1 - k • v 



(5) 



a(*o) ll/2 



a(t) 



b(i ) 



which consists of three factors. The first factor is the special relativistic Doppler shift 
while the second and the third factors are general relativistic corrections. They are 
rather complicated functions of positions and velocities of observer, spacecraft, and 
light-ray-deflecting masses. Exact expressions for these functions can be found in 



7 Post-Newtonian Theory for Doppler Measurements of Binary Stars 

The determination of velocities of stars from precise Doppler measurements is de- 
scribed in We apply successive Lorentz transformations and the relativistic equa- 
tion of light propagation to establish the exact treatment of the Doppler effect in 
binary systems both in Special and General Relativity theories. As a result, the 
Doppler shift is a sum of (1) linear in c _1 terms, which include the ordinary Doppler 
effect and its variation due to the secular radial acceleration of the binary with respect 
to observer; (2) terms proportional to c -2 , which include the contributions from the 
quadratic Doppler effect due to the relative motion of binary star with respect to the 
Solar system, motion of the particle emitting light, orbital motion of the star around 
the binary's barycenter, diurnal rotational motion of observer, and orbital motion of 
the Earth; and (3) terms proportional to c~ 2 , which include the contributions from 
redshifts due to gravitational fields of the star, star's companion, Galaxy, Solar system, 
and the Earth. We briefly discuss in the paper feasibility of practical implemen- 
tation of these theoretical results, which crucially depends on further progress in the 
technique of precision Doppler measurements. 

The author thanks G. Neugabauer for constant support and G. Schafer for useful discussions. 
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